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Abstract. From the analysis of Dopplergrams in the K i 7699 A and Na i 5890 A spectral 
lines observed with the Magneto-Optical filter at Two Heights (MOTH) experiment during 
the austral summer in 2002-03 we find upward traveling waves in magnetic regions. Our 
analysis shows that the dispersion relation of these waves strongly depends on whether the 
wave is detected in the low-beta or high-beta regime. Moreover, the observed dispersion 
relation does not show the expected decrease of the acoustic cut-off frequency for the field 
guided slow magnetic wave. Instead, we detected an increase of the travel times below the 
acoustic cut-off frequency and at the same time a decrease of the travel time above it. To 
study the formation height of the spectral lines employed by MOTH in greater detail we are 
currently in the process of employing 3D MHD simulations carried out with CO^BOLD to 
perform NUTE spectral synthesis. 
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1. Introduction 

The physics of the chromosphere and corona is 
most important for an understanding of the for- 
mation and variability of the solar UV/EUV ra- 
diation. However, the heating mechanism that 
causes the rising temperature profile of the up- 
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per chromosphere is not fully understood. In 
particular, the reason why the Sun shows a ris- 
ing temperature profile over all phases of the 
solar cycle is still an open question. 

A number of potential heating mecha- 
nisms of the solar chromosphere and coro- 
na are e.g. low-freque ncy magnetoacoustic 
gravity (MAG) waves dJefferies et all l2006i: 
iMcIntosh & JefferiesLl2006ir the Farley-Bune- 
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Fig. 1. Observed dispersion relations from the anal- 
ysis of the Doppler signal observed with MOTH. 
The thin lines denote the phase travel time and 
the thick lines the group travel times. The differ- 
ent linestyles denote whether the travel time was 
derived from a magnetic or non-magnetic region, 
i.e. the dotted line corresponds to areas on the solar 
disk where a < c, the dashed lines where a ~ c, 
and the solid line where a > c. For details see 
iHaberreiter & Finsterld fcOlOh . 



man instability jFontenla et al.L l2008l). coronal 



heatin g from Alfv en waves |Cranmer et al 



(2007), nanoflares iPatsourakos & Klimchuk 
(|2009|) . and heating processe s of the corona are 
rooted in the chromosphere dPe Pontieu et al.L 
l2009h . 

The ultimate aim of this work is to study 
the energy provided by MAG waves. Theory 
predicts that different types of MAG waves ex- 
ist in t he solar atmo sphere jSchunker & Callvl 
120061; ICallvl l2007l) . They can couple at the 
height where the magnetic pressure equals the 
gas pressure, exchange energy, and thereby 
convert into a different wave type. Moreover, 
due to partial mode conversion, theory predicts 
that more than one wave type is likely to ex- 
ist on both sides of the conversion layer. These 
processes are already well understood from a 
theoretical point of view, however mode con- 
version has not yet been unambiguously iden- 
tified in observations. 

The unambiguous identification of the 
different wave types from the observations is 
however a challenging task. To clearly disen- 
tangle the various wave types it is crucial to 



determine the observational heights of the in- 
strument with respect to the ma gnetic canopy, 
as ha s been pointed out bv IBogdan et alj 
(l2003l) . From the Mag neto-Optical Filter at 
Two Heights experiment dFinsterle et al.Ll2004l 
MOTH) and the Michelson Doppler Imager 
observations dScherrerl Il998l MDI) onboard 
the NASA/ESA So lar and Heliospheric 
Obser vatory (SOHO), IHaberreiter & Finsterl j 
d2010l) determined the group and phase travel 
times of the observed MAG waves depending 
on the a - c-layer, i.e. the layer where the 
sound speed equals the Alfven speed, with 
respect to the observational height of the Ki 
7699 A and Nai 5890 A lines. Fig.Q gives 
the dispersion relation for the cases where 
a > c, a ~ c, and a < c ; for details see 
IHaberreiter & Finsterld d2010l) . 

From this analysis it is clear that the 
conversion layer changes the type of wave 
that is observed. The main question that 
needs to be answered is to what extent 
is the wave converted from a slow acous- 
tic m ode to a fast magnetic mode dCallvL 
l2007h . The observed dispersion relation has 
the potential to shed light on this matter. 
However, the dispersion relations derived from 
the observations can only be fully under- 
stood through the analysis of wave propa- 
gation from numerical simulations. We will 
use 3-dimensional (3D) MHD sim ulations car- 
ried out with CO ^BOLD (Frevtaa etaL, 2002'; 
Wedemever et ali ,2004.; Schaffenberger et al.i . 



2006HSteiner et alil2007l)~ 

In this paper we report on continuing ef- 
forts to provide synthetic data that might allow 
us to understand the observed dispersion rela- 
tions of these waves. 

2. Ongoing spectral synthesis from 
3D IVII-ID simulations 

Fig.|2|shows the temperature profiles of a snap- 
shot of a vertical slice through the simulation 
carried out with CO^BOLD. We are currently 
in the process of calculating the NLTE level 
populations of hydrogen, sodium, and potas- 
sium, and the most important elements to syn- 
thesize the solar line profiles for each of the 
grid points and snapshots. This will allow us 
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Fig. 2. Panel (a): Shown are the temperature structures of a snapshot of the simulation carried out with 
CO'^BOLD. Panel (b): Shown is the modulation of the K i 7699 A line calculated with COSI based on 
the Doppler signal of a vertically propagating wave with v=14mHz in a ID atmosphere structure. The 
dashed line indicates th e narrow band MOF filters employed with the MOTH instrument; for details see 
iHaberreiter et alj j2007r) . 



to study in detail the formation height of these 
lines, which is essential for an understanding 
of the phase lag between the different observa- 
tional layers in the solar atmosphere. 



3. Results from 1D model 

For testing reasons we have akeady modeled 
the K I 7699 A line profile being modulated 
with the Doppler velocity signal correspond- 
ing to waves with various frequencies. The 
line profiles have been ca lculated with the 
code for Solar Irradiance dHaberreiter et al.L 
120081 COSI). Fig.|2t) shows the line profiles for 
waves with 14mHz. It is clear that the spectral 
line does not show a pure Doppler shift, but 
rather a change in the overall line profile, i.e. 
the C-shape of the line. The reason is that the 
core of the line is formed higher in the solar at- 
mosphere than the wings of the line. Therefore, 
the line center and the wings of the line are 
not modulated simultaneously by the Doppler 
signal but experience a time lag. If the spatial 
wavelength of the acoustic wave is of the order 
of the extent of the contribution function of the 
spectral line, then the line profile is modulated, 
as shown in Fig.|2j5. This poses a challenge 
for the exact determination of the Doppler sig- 
nal of high frequency waves from this kind of 
observ ations; for details see IHaberreiter et alj 
(l2007h . 



Detailed modeli ng of the wa v e exci tation 



as carried out bv iSteiner et a l. (2007) and 



iNutto et aP (lin press, see this volume.) will al- 
low us to disentangle the effect of the wave- 
length of the acoustic wave on the line profile. 



4. Synergy of 1 D and 3D models 

Semi-empirical ID models are very success- 
ful for modeling the solar spectrum from 
the UV/EUV to the IR (Fontenlaet al., 2009; 
IHaberreiter et aU 12001 IShapiro et all l2010l) . 
Of course, ID models do not represent any dy- 
namics of the solar atmosphere, but are a tem- 
poral and spatial mean of the solar atmosphere. 
However, as they can reproduce a broad wave- 
length range and also detailed line profiles, 
they are a powerful means to study realistic ra- 
diative losses of different layers in the solar at- 
mosphere, in particular the chromosphere, for 
the quiet Sun as well as active regions. 

In the near future we plan to synthesize de- 
tailed line profiles as well as the spectrum over 
broad wavelength ranges based on 3D MHD 
simulations. Hopefully, ID and 3D models can 
provide complementary insight to achieve this 
goal. For example, detailed radiative losses cal- 
culated from ID models might prove valuable 
for a comparison with the values used in 3D 
MHD simulations. 
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5. Discussion and Conclusions 

From the analysis of Dopplergrams in the 
K I 7699 A and Na i 5890 A spectral lines 
observed with the MOTH experiment dur- 
ing the austral summer in 2002-03 we find 
upward traveling waves in magnetic regions. 
Our analysis shows that the dispersion relation 
of these waves strongly depends on whether 
the wave is detected in the low-/3 or high-yS 
regime. Moreover, the observed dispersion re- 
lation does not show the expected decrease 
of the acoustic cut-ofF frequency for the field 
guided slow magnetic wave. Instead, we de- 
tect an increase of the travel times below the 
acoustic cut-ofF frequency and at the same time 
a decrease of the travel time above it. To study 
the formation height of the spectral lines em- 
ployed by MOTH in greater detail, we use 3D 
MHD simulations carried out with CO^BOLD 
to carry out NLTE spectral synthesis. 

Total and spectral solar irradiance changes 
are understood to be linked to the Sun 's mag- 
netic cycle; see e.g. lHaberreitej (12010 ). With 
increasing magnetic activity more radiation is 
emitted, in particular in the UV and EUV. 
However, the physical processes that pro- 
vide the enhanced heating in magnetic regions 
are not yet fully understood. Waves leaking 
through magnetic regions are considered one 
of the drivers. A detailed spectral synthesis 
based on different scenarios of magnetic field 
strengths representing different phases during 
the solar cycle might help to understand the 
contribution of MAG waves in the heating pro- 
cess. 
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